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Abstract 

A theory of anyon excitons consisting of a valence hole and three quasielec- 
trons with electric charges (— e/3) is presented. A full symmetry classification 
of the k = states is given, where k is the exciton momentum. The energy 
levels of these states are expressed by quadratures of confluent hypergeomet- 
ric functions. It is shown that the angular momentum L of the exciton ground 
state depends on the distance between electron and hole confinement planes 
and takes the values L = 3n, where n is an integer. With increasing k the elec- 
tron density shows a spectacular splitting on bundles. At first a single anyon 
splits off of the two-anyon core, and finally all anyons become separated. 
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I. INTRODUCTION 

Recent progress in the intrinsic interband spectroscopy of the Fractional Quantum 
Hall Effect (FQHE) || made a demand for developing a theory of optical phenomena ac- 
counting in an appropriate way the basic properties of FQHE phases. Such a theory has to 
reflect two essential facts: i) an exciton nature of the low-energy states in a strong magnetic 
field, and ii) the presence of an Incompressible Quantum Liquid (IQL) || underlying the 
FQHE. Therefore, the theory should describe an exciton against a background of an IQL. 
There are two approaches to this problem. The first one is based on exact numerical treat- 
ment of few-electron systems in a spherical geometry |7HT2f . The second one is based on the 



anyon exciton (AE) model |13| which considers an exciton as a neutral entity consisting of 
a valence hole and several quasielectrons (QE's) carrying fractional negative charges. These 
both approaches have a restricted applicability depending on the value of the parameter h/l, 
here h is the separation between electron and hole confinement planes and I = (ch/eH) 1 ^ 2 
is the magnetic length. Reliable calculations in a spherical geometry are accessible only for 
h/l < 1.5. On the contrary, the AE approach fails when h/l is small, but is exact in the 
opposite limit h 3> I. Indeed, by electrostatic arguments the size of an AE ground state is 
about h, and the model is applicable when the mean distance between anyons is large as 
compared to the anyon size which is about l; in the same limit the perturbation exerted by a 
hole is small. There is our believe that both approaches will result in a similar picture in the 
intermediate region h/l ~ 1 and, therefore, taken together they will provide a description 
of excitons in the whole region of the h/l values. Experimental data on excitons are not 
avalable for the region h > I yet, and therefore a comparison with the experiment can not 
be performed at present. For this reason we abstain from discussing the involved problem of 



the transition probabilities |L|] and concentrate on the energy spectra and electron density 



distribution. The relation between the results of the AE approach and of computations for 
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few-particle systems is discussed in the last Section. 

The oversimplified model of an AE consisting of a hole and two semions (anyons having 
charges (— e/2)) has been developed by us previously ||13|| . In that paper a multiple-branch 



spectrum of an AE has been established as a basic signature of the charge fractionalization 
in IQL's. Recent data [I3|,r5| have revealed the appearance of this structure also in few- 



particle computations. It is remarkable that several new exciton branches appear at very 
moderate values of h/l % 1. In what follows we present for the first time a theory for a 
realistic case of three QE's which corresponds to a v = 1/3 IQL. 

II. THEORY 

Let us consider an exciton consisting of a valence hole with a charge (+e) and three QE's 
with electrical charges (— e/3) and statistical charges a = —1/3 (a = for bosons and a = 1 
for fermions). In the strong H limit, when the Coulomb energy Ec = e 2 /el <C hu c , where 
u c is the cyclotron frequency, it is convenient to use dimentionless variables scaled in units 
Ec, I and e. Instead of hole, r^, and anyon, r*, coordinates it is convenient to introduce the 
following two-dimensional coordinates: 

1 1 3 1 3 

R = ( T h + o X, ' P = o X, r i ~ r h > T ij = T i ~ r j> *j = 12 > 23 > aIld 31 ' W 

2 d i=l d i=l 

Complex coordinates Zji = Xji + iyji, as well as r^, are not independent. Indeed: 

ri2 + r 23 + r 3 i = . (2) 

Since AE is a neutral entity it possesses an in-plane momentum k [|l(|. As a basis for 
the exciton wave functions in the Halperin's quasi wave function representation [17,18] the 
functions 

#L,k = exp{*kR + i( Px Y - p y X) - (p - d) 2 /4}P L (. ..z jt . . .) I](^) a exp{-|^ 7 | 2 /36} (3) 

ji 

may be chosen. Here d = z x k, and (— d) is the exciton dipole moment. Exponential factors 
in Eq. |3] ensure the translational magnetic symmetry and belonging the function to the 



ground Landau level both for anyons and a hole. A function Px(- ■■%••■) is a homogeneous 
polynomial of the degree L which is symmetric in all coordinates Zi. The function ^ 
describes an exciton with the momentum k and the internal angular momentum (-L). 

For a three anyon problem polynomials Pl depend on three variables zu, Z23, and z%i. 
If to take into account the constraint Z12 + z 2 % + £31 — 0, one can show that the total number 
of linearly independent symmetric polynomials equals [L/6] + 1 for even L and [(L — 3)/6] + 1 
for odd L. Square brackets designate the integer part. Odd L polynomials start with L = 3 
because of Eq. |[ Even-L polynomials Pl may be classified as: 

p _ =L-4M=2M=2M , -L-AM -2M -2M , -L-4M -2M -2M / a\ 

r L,M — *12 ^23 ^31 ^23 ^31 ^12 "r ^31 ^12 ^23 > VV 

where M = 0, 1, ... [L/6]. The only symmetric polynomial with L = 3 is a Vandermonde 
determinant in the variables Zjf. 

V = {z~i2 - ^23) (^23 - %) (231 - ^12) • (5) 

As distinct from the polynomials Pl,m with even L, which are symmetric both in the bosonic 
permutations z% <-> ^2 and permutations 212 <-> ^23, etc., the polynomial V is symmetric in 
the permutations of the first and antisymmetric in the permutations of the second type. All 
independent odd-L polynomials can be represented in the form: 

Pl,m = VP L - 3tM , Ps, = V. (6) 
To our best knowledge, in the previous studies only the even-L polynomials have been taken 



into account [fL9|| . When choosing polynomials Pl,m we have not imposed the hard core 
constraint and defer the discussion of the related properties to what follows. 

In the basis of polynomials Pl,m the set of quantum numbers consists of k, L and M. The 
only disadvantage of the functions \& 'l,mm is that the scalar products < LMk|LM'k >^ 
for M 7^ M'. As a result, the matrix B of these scalar products which is diagonal in k 
and L has a block- diagonal form; the block sizes increase with L. The Schrodinger equation 
is Hx = tP>x- All integrations in scalar products and matrix elements of H should be 



performed over four variables r^, ri, r2, and r3. It is convenient to choose the new variables 
R, p and three and take into account the constraint of Eq. ^| by the usual 5-function 
transformation. It adds a new variable f, but all calculation become symmetric in the 
anyon variables. As an example we show here the transformation of a matrix element with 
polynomials Pj_,,m substituted by monomials z™| ^23^31- After the 5-function transformation 
the matrix element takes the form: 

f df 

< n 1 n2n 3 k|n 1 n 2 n 3 k >= J ^^p-^nm^ 

(f)AC 2n ,(f)A^ 3n ,(f) , (7) 



where 



Ad" (f) = 27Ti' n_n '' -^( maX { ?7 '' n '} + Q + 1) 2("+™')/2+a ^{n+n')+2(a+l) ^\n~n'\/2 

\n — n'\\ 

x exp{i(pf (n' — n)} $(max{n, n'} + a + 1, \n — n'\ + 1 ; —t) . (8) 

Here t = 9f 2 /2, iff is the azimuth of f, and $ is a confluent hypergeometric function. 

Quasi wave functions of Eq. |3] differ from bosonic functions only by a statistical factor 
Ylji(zji )~ a . It is of most importance for small values of |zy| ~ I. For the same distances the 
deviation of the anyon-anyon interaction law from a Coulomb law [2C] also should be taken 
into account and results in comparable contributions. In what follows we restrict ourselves 
with a Coulomb interaction in the Hamiltonian H = V aa + V a h , where V aa = (1/9) Y,ij T Tl 
and V a h = — (1/3) J2i \ v m + hz\~ x , and omit the statistical factor in ^l,m,\l- This last 
approximation strongly simplifies all calculations. Indeed, with a = all functions $, 
Eq. turn into polynomials if one performs the Kummer transformation, $(/3, 7; t) = 
e*$(7 — /3, 7; — t) and takes into account the fact that 7 — (5 are non-positive integers 
for all functions (^). As a result, integrals (J8|) can be calculated exactly. The matrix 
elements of V aa differ from those of Eq. [7| by changing in one of the Ai^ n i functions a 
for (a — 1/2). This latter function can not be reduced to a polynomial, hence, one-fold 
integration over / should be performed numerically. The most complicated are the matrix 
elements of V a h ■ To simplify them it is convenient to perform a Fourier transformation 
V a h( r ) — I rfqV^q) exp(iqr)/(27r) 2 . In this representation the integrand contains again 



three factors M°t-n' • O ne of them has f as an argument, while two others have (f ± q/3) . 
For a = the integrations over / and over the angle between f and q may be performed. 
We choose ^l,m,u in such a way that the matrix elements of H are real, and the diagonal 
elements of B equal unity. The final expression for the matrix element of V a h is: 

POD 

< LMk\V ah \L'M'k >= - / exp(-3g 2 /2 - qh)J\ L -i/\{kq)Q LMiI/M ,{q)dq , (9) 

Jo 

where J\L-u\{kq) are Bessel functions, and the functions Qlm,l'm'(q), rea l an d symmetric 
in the LM indices, are polynomials in q. If k = 0, the integral equals zero for L ^ V 
as it follows from the angular momentum conservation law. Therefore, for k = the 
total Hamiltonian has a block-diagonal form. The lower polynomials are of a simple form: 
Qoo.oo = 1, Q2o,oo = <? 2 /2, ^20,20 = I - q 2 + q 4 /4:. Since the blocks are of the 1 x 1 
size for the even-L polynomials with L < 6 and the odd-L polynomials with L < 9, the 
corresponding wave functions do not depend on h, and the energies may be expressed in terms 
of quadratures from transcendental and elementary functions. If L is even and 6 < L < 12, 
the coefficients of 2 x 2 secular equations may be expressed in the same terms, etc. With 
increasing h the anyon-hole attraction becomes smaller, the effect being stronger the less is 
L. The block-diagonal structure of H allows the energy level intersections for k = 0, for 
k these level crossings turn into anticrossings. 

In Fig. the energy spectrum e(k) is shown for two values of h. The following regu- 
larities are distinctly seen. With increasing h, the levels with higher L values draw closer 
to the spectrum bottom. The branches with small L values have a larger dispersion. The 
anticrossings are wider the larger the k values and the lesser the differences \L — L'\ are. It 
is remarkable that the branches with a negative dispersion near k = exist; it may appear 
even in the ground state, Fig. |l|a. For a two-semion problem the dispersion near k = is 



always positive [13 



The most interesting property of AE's which follows from our calculations is the distri- 
bution of the anyon density Di(r, k) around the hole 
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A(r,k) = — £ cos[(L-L')0] 

Z7F LM,L'M> 

x dqq exp(-3g 2 /2) J, L _ L ,|(g|r - d|) Qlm,L'M'(q) Xi*M'( k )xLtf(k) , (10) 

JO 

where Z numerates spectrum branches, and 9 is the angle between the vectors d — r and d, 
d is the center-of-mass of the anyon density. 



III. ZERO MOMENTUM ANYON EXCITONS 

For k = the density distribution D(r, 0) is shown in Fig. ||] for the states with L < 6. 
It has been mentioned above that Dl do not depend on h for L < 5. The state L = 3 
is remarkable in some sense. It is the first to show a density minimum at r = 0, and it 
generates all even-L polynomials, Eq. |6|. For L = 6 there are two eigenfunctions; they 
depend on h. In Fig. |] they are shown for h = 0; the lower energy component is drawn by 
a solid line. With increasing L the functions Dl^m become broader, and as a result some 
of the states with large L values draw closer to the spectrum bottom when h is growing, 
Fig. |T[ However, it is seen from Fig. |3] that only the states with L = 3n, n > are integers, 
reach the spectrum bottom. The bottom states described by even- and odd-L polynomials 



alternate. The periodicity with the period AL = 3 resembles the Laughlin's result |2T 
for a three-electron system. We attribute this periodicity in L to the quantization rule for 
identical particles, since the rotation angle between two exchange points (which enters into 
the semiclassical quantization rule |22j) equals 2tt/3. 



The functions with L = 6M, M > 1 are integers, depend on h strongly. For L = 6 
the electron density distribution D(r, 0) changes rapidly in the vicinity of h — 2, Fig. |j. 
For h < 2 the low energy component of the L = 6 state is close to \l/6,o (fc> r h = the 
overlap is 0.96), while for h > 2 it approaches the hard core state \I/ 6 ! (for h = 3 the 
overlap is 0.97). The function \l/6,i is the first in the series of the even-L polynomials 
Pgm,m = (%2^23^32) 2M , M > 1, forming ground states for the large h values. The density 
-Dl,m(^, 0), L = 6M, has a single maximum for each of these states. The maximum of 
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the density Dl,m{t-, 0) found after averaging over is determined by a simple equation 
tl = V%L', the maximum of Dl,m{t,^) is pretty close to this point. In the vicinity of 
the maximum of -Dl,m(?", 0) the three particle correlation function reaches the maximum 
for an equilateral triangle. It is just the result which is expected in the L > 1 limit from 
electrostatic arguments. The value of r?-, found both as a mean value in a ^gm,m state and 



by the semiclassical approach in terms of the distances between anyon guiding centers |f22 
satisfies the equation L = 6M = r\ 2 /& — 2. 

The density distribution for ^3^(r) shows a single maximum, Fig. 2. In the vicinity 
of it the three-particle correlation function reaches the maximum for an equilateral triangle 
configuration of anyons with r\ 2 w 18. In this respect the functions Pg,i and P^^ show a 
similar behavior. A strong difference between them is reflected in the anyon two-particle 
correlation functions W(r): 

= ^(1 + ^(r/6) 4 )exp(-r 2 /12), 1^,0 = ^(1 + \{r/2f) exp(-r 2 /12). (11) 

They are shown in Fig. |5| . The function We,i has a hard core behavior, while W^p reaches 
its absolute maximum at r = and has the second maximum which is only by a factor 
0.97 lower than the main one. Of course, all bottom state odd-L polynomials with L > 9 
show a hard core behavior. It follows from the above data that with increasing h the charge 
fractionalization results in the broadening of the anyon density distribution in the low-energy 
k — states. The dependence of the exciton binding energy on h is close to the Coulomb 
law, but the numerator is considerably less than unity. In Fig. the energy of an ordinary 
k = magnetoexciton, EmeW, is shown for comparison. Because of the existence of a 
hidden symmetry, Sme{J^) gives an exact result in the h — > limit. The AE model becomes 
exact in the opposite limit, h — > oo. The both curves should be matched somewhere in the 
intermediate region, h ~ 1. 
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IV. MOMENTUM DEPENDENCE OF THE EXCITON CHARGE DENSITY 

Since the exciton dipole moment d differs from k only by the factor el 2 and the rotation 
by 7r/2, Sect. [TT], one can expect that with increasing k the electron density splits into 
bundles, their charges being multiples of e/3. The splitting of the electron shell into two 
well separated quasiparticles has been observed previously for a two-semion exciton ||13|| . For 
a three-anyon exciton the pattern are much more impressive. They are shown in Fig. |B| for 
h = 3 when the criterion of the large electron-hole separation is fulfilled. The distribution 
which is cylindrically symmetric for k = transforms with increasing k into a distribution 
with a single split off anyon, k = 2 and 3. Two anyons constituting the exciton core show a 
slight but distinct splitting in a perpendicular direction. This core changes its shape with k, 
but remains stable for a long. Finally, for rather large k values, it splits in the d direction as 
it is seen in the last figure, k = 6. In all figures the center-of-mass of the electron density is 
at the point x = k, and the asymmetric distribution of the density arises due to the odd-L 
polynomials. In the absence of them the distribution had to be symmetric with respect to 
the point d. 

The well outlined profiles of the electron density seen in Fig. |^ may be somewhat 
smeared by the oscillatory screening inherent in IQL's J53|. Nevertheless, the basic pat- 
tern of the charge separation in an exciton should strongly influence the k dependence of 
the magnetoroton-assisted recombination processes since charge density excitations are left 
in a crystal afterwards. 

V. DISCUSSION 

The above theory is based on a Coulomb interaction between all the particles and ne- 
glecting the statistical factor T[ji(zji) a m Eq. |3|. For small L values the theory results in 
dense states which can not be described by the AE model, and for larger L values in more 
loose states to which the AE model should be applicable. The L3 state or, more probable, 



9 



the Lq state are the first candidates. If one takes into account the statistical factor, the 
effective repulsion of anyons should increase. Therefore, we expect that the critical value 
of h, h CT , when the loose states reach the spectrum bottom must decrease. According to 
Fig. [I h CI pd 2 in the AE model, while the computations in the spherical geometry result in 
h CT ~ 1 for v = 1/3 |H|jr5 |. These data are in a reasonable agreement. The more intriguing 



question is whether the realistic anyon form-factors and quasipotentials may result in some 
new types of the k = ground states, e.g., the two-scale states like the low-energy compo- 
nent of the L = 6 state in Fig. |], h = 1. It has been proposed by Apalkov and Rashba ||, 
that the two-scale exciton observed by them in computations for a v — 1/3 IQL consisted 
of a two-anyon core and a split off anyon (it resembles the exciton shown in our Fig. || 
k = 2, if described in the (L,k) rather than k representation). Such a configuration may 
be energetically favorable for a k = exciton if the two-QE repulsion at small distances is 
suppressed well below its Coulomb value. This hypothesis is in agreement with the data of 
Beran and Morf ||20|| . Therefore, the configurations of anyon shells of AE's with moderate 
values of h are expected to be sensitive to the QE interaction law, and calculations based 
on realistic quasipotentials are desirable. 
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FIGURES 

FIG. 1. Anyon exciton dispersion law e(k) for two values of h. Mention the opposite signs of 
the ground state dispersion for both curves. Anticrossings become tiny with increasing h. Numbers 
show L values. For more detail see text. 

FIG. 2. Axisymmetric electron density distributions Dl(v, 0) for the states with L < 6. Two 
L = 6 states are shown for h = 0; the lowest state is shown by a solid line. Numbers show L values. 

FIG. 3. The energy e(0) and the electron density -Dl(0, 0) at the point where the hole 
resides plotted vs h for the ground state of an exciton with k = 0. The ground state energy of 
an anyon exciton is shown by a solid line; the dots on it show the positions of the intersections 
between the energy levels with different L-values. For comparison the energy of a conventional 
magnetoexciton eme with k = is shown by a dashed line. Numbers near the Dl(0, 0) curve show 
the L values. Only the states with L = 3n reach the spectrum bottom (as an exclusion the state 
L = 2 appears as a bottom state in an extremely narrow region of the h values). 

FIG. 4. Electron density distribution D(r, 0) in the L = 6 states with k = for three values of 
h. The density in the lower energy state is shown by solid a line. Consecutive numbers, I, of the 
energy levels are shown near the curves. 

FIG. 5. Normalized anyon pair correlation function W(r) for the states *6,i and ^3,0; k = 0. 

FIG. 6. Electron density distribution in an anyon exciton for different values of k. A hole is in 
the origin, the x axis is chosen in the d direction. The center of the electron density distribution 
is at x = k. 
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